tATlONAL ADVISORY COMMlTlEE 
. FOR KERONAUTICS 

iUN 20=924 TECHNICAL NO^ES 

MAILED NATIONAL ADVISORY OOlllflTTEE K)R AERONAUTICS. 



No. 196 



REMARKS ON THE PRESSURE .DISTRIBUTION OVER THE SURFACE OF AN 
ELLIPSOID, MOVING TRANSLATIONALLY THROUGS A PERFECT FLUID. 

By Max M. Munk, 



June, 1924. 



NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS. 



TECHNICAL NOTE NO. 196. 

REMARKS ON THE PRESSURE DIST?."CETTTI-ON OTIHR THE SURFACE OF AN 
ELLIPSOID, MOVING TRANSLATIOil-ALLY OiiROUGH A PERFECT FLUID. 

By Max M. Mank* 

Summary 

This note, prepared for the National Advisory Committee for 
Aeronautics, contains a discussion of the pressure distribution 
over ellipsoids when in translatory motion through a perfect fluid. 
An easy and convenient way to determine the magnitude of the veloc- 
ity and of "the pressure at "each point of the surface of an ellip- 
soid of rotation is described. 

The Imowledge of sudi pressure distribution is of great prac- 
tical value for the airship designer. The pressure distribution 
over Ihe nose of an airship hull is Imotm to be in such good agree- 
ment with the theoretical distribution as to permit basing the com- 
putation of the nose stiffening structure on the theoretical dis- . 
tribution of pressure. 

References 

1. Horace Lamb - "Hydrodynamics," Chapter V. 

2» Max M. Munk - "The Aerodynamic Forces on Airship Hulls," 

N.A.O.A. Technical Report No. 184. 

3. R. Jones & D. H. Williams - "The Distribution of Pressure 

Over the Surface of Airship Model N 721," 
Reports and Memoranda No- 600. 



N.A. C-A. Teclmical Note No- 196. " 2 

■ Experiments have shown that the knoiiTiedge of the pressure dis- 
trilmtion over the surface of ellipsoids, moving trahslationally 
through a perfect fluids is often of considerable practical inter- 
est. This pressure distribution is of a simple description; it 
can easily and quickly be determined by anal3rtical methods. To 
the best of my knovrledge this has never been brought out clearly 
in amy publication. The mathematical theory of the flow created 
by an ellipsoid is given by H. Lamb in his "Treatise on Hydrody- 
namics,'^ Chapter V. It requires considerable mathematical train- 
ing to grasp the full meaning of the results as given by Lamb. 
E. Gallop (Ref . 3) has given some comments on the nature of the 
resulting distribution of the velocity and pressure, for the spec- 
ial case of an ellipsoid of revolution. A part of this holds true 
for all ellipsoids, including those with three different principal 
axes. Mr. Gallop does not, however, for the special case of sphe- 
roids make 12ie distribution of the pressure sufficiently plain for 
immediate computation or for the practical application of this 
interesting analysis. 

The knowledge of one simple lemma on the potential flow around 
ellipsoids, implicitly contained in Lamb's result (Third Edition, 
Equation (114 (8) ), ) is sufficient for the deduction of all the 
following theorems and for the determination of th'e pressure dis- 
tribution. This lemma is^ 

If an ellipsoid is moving with uniform velocity parallel to 
one of its principal axes, say parallel to the x-axis, the velocity 
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potential at any point of the surface can be written in -ttie form 



where A* is constant for a given flovf and a given ellipsoid. 
This theorem is the key to all the relations referring to the d is- 
tribution of velocity and of presEJure. 

If the velocity of flow is not parallel to a principal axis, 
but has components in the direction of each of them, the resulting 
flow is the superposition of three flows analogous to the one just 
considered. Hence, at all points of the surface, the .potential is 
a linear function of the Cartesian coordinates x, y, and z again, 
and can be written in the form 



TAiere the coordinate axes are chosen to coincide with the axes of 
the ellipsoid. Hence the curves of equal potential ^ are situated 
on parallel planes. 

KTow, suppose first the ellipsoid to be at rest and the fluid 

» 

to be moving relative to it, as in a wind tujonel or as with an 
airship moored in a gale. The change from the ellipsoid moving 
through the fluid otherwise at rest to the fluid passing by the 
stationary ellipsoid does not affect the validity of Equation (s) 
except giving the constants A», B«, and 0* other values, say 
A", B", and 0". In the latter case (the body at rest) the velocity 
of the fluid at all points of the surface is parallel to the sur- 
face. Consider first the elements of surface ^containing a line 



^ = A'x 



(1) 



$ = A»x + B»y + O^z 



(2) 
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element at right angles to the planes of constant potential, i^e* 
at the points of the ellipsoid where the plane A"x + B"y + G"z = 0 
meets the surface. It is apparent from (3) that at all these points 
the velocity has the coniponents A"', B", and 0". This is e^ridently 
the mx3.imim velocity. 

At all other points of the ellipsoid the elements of surface 
are inclined towards the direction of raaxiraam velocity, say "by the 
angle e . Then the elements of distance on tbe surface, As, be- 
tween curves of eqiial notential are increased in the ratio — \ — , 

cos € 

when compared with the actual distances hetween the planes of equal 
potential. Accordingly, the velocity, being equal to is de- 

creased inversely, its mgnitude is A" cos c It will be noted 
in particular that the velocity is equal at surface elements whi ch 
are inclined by the same angle c . It is equal to the projection 
of the maxinram velocity at right angles to tbe surface element- 
Hence the velocity cannot exceed the one rightly denoted by "maxi- 
mum velocity," having the components A", B", and 0". 

Returning to the case when the direction of flow is parallel 
to a principal axis, it can be shown that the maximum velocity A" 
sta.nds in a very simple relation to the kinetic energy of the flow, 
and hence to the apparent additional mass of the ellipsoid- We 
have now to suppose the fluid to be at rest and the ellipsoid to 
move, say with the velocity U, parallel to a principal axis, 
e»g», the x-axis« The kinetic energy of the flow set up is equal 
- g'P/^^dS, i.e. the volume of fluid displaced by an 
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element of tlie surface per unit of time, multiplied by the poten- 

p 

tial at the point of dispDaoeiaent and by g where p denotes the 
density of the fluid. Sfow, the volume displaced by a surface ele- 
ment per unit time is equal to -Hie projection of this element 

perpendicular to the direction of x, multiplied by the velocity 

P 

tJ. The potential being A'x the integrand becomes A- Ux g dy dz. 

/ xdy dz is the volume of the ellipsoid, hence the integral gives 

Volume A*"!! ^. This is the kinetic energy, usually expressed by 
— ^£ 

Volume , where kj denotes the factor of apparent mass* 

It follows that 




A", referring to the case when the ellipsoid is stationary, is 
connected with A*' by the equation 

A" = A' + Ui 

as the latter flow results from the former by the superposition of 
the constant velocity U» Hence it appears that 

A!L=k, + i = A. 

A is the maximum velocity corresponding to a flow having unit 
velocity along the x-axis. It is a constant for a given ellipsoid. 

It equals the sum of 1 and of the factor of apparent mass 
as is confirmed for two special cases, where the factor A is well 
known. With a sphere, the maximum velocity is 1»5 times the veloc- 
ity of flow, and the additional apparent mass is one-half the mass 
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of the displaced fluid. With a circalar cylinder, moving at right 
angles to its axis, the aiaxiimirf. velociby is tiivics the velocity of 
flow, and the apparent additional mass is 3-:n:9l to the mass of the 
displaced fluid. 

A" B" 0" 

The iratios A = |j— , B = y-, ^ = %~ are independent of the 
velocities U, V, and W and hence only depend on the ratio of 
the three semi-axes of the ellipsoid' a, b, and c. Lamb gives the 
method to compute them. Oonpite first the integral 



* dx 
a = a b c / ^ 

o 



(a= + x) y (a^ + x) (b^ -i- x). (c^ + x) 



and the analogous integrals 3 and 7 for the axes b and o. The 
factors of apptarent fijass are then 

There are no tables for A, B, C, or for kj, kg, kg, pub- 
lished yet. The integrals for a etc, can be numerically evalu- 
ated in each case, and I will assume at present that A, B, and 0, 
are therefore known- For the special case b = c, that is, for 
ellipsoids of revolution, k^ and kg have been computed for a ser- 
ies of elongation ratios,' and are reprinted in a small table in 

Ref . 3. They are connected by the relation k, = . The 

2 kg 

determination of the velocity at any point is thus reduced to a 
simple geometric problem. The maximum velocity, whose components 
are AU, BV, and CVSi\ has to be projected onto the plane tangent 
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to the ellipsoid at the point considered, i.e. it has to "be multi- 
plied by the cosine of the angle "betinfeen the normls to the oux- 
face at this point and at the point wisve tie velocity is a maximum- 

lik the most interesting case of an ellipsoid of revolution, 
this can' he done analytically in a very convenient way- The for- 
mula is most easily arrived .at by the application of elementary 
vector analysis. First compute the component of the maximm veloc- 
ity in a direction normal to the surface at a given point. The 
longitudinal component of the maximam velocity is AU, and the 
lateral component of the maximal velocity. is BV. 

Let the angle between the normal and the longitudinal axis be 
f] , and let the dihedral angle between the plane containing this 
axis and the line of velocity of flow and the plane containing the 
axis and the point in question be p. Then cos ri is the longi- 
tudinal component and s±nr\ cosP the lateral component of the 
no"rmal of unit length. Hence the component of the maximum velocity 
in a direction perpendicular to the element of surface is 

= (1 + ki) U COST] + (1 + ka) V sinT) cos 

Let Vj^ denote the component parallel to the surface element. 

Then 

+ Vg^ = '^'max* hence 
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This is tile desired fornrala for the velocity of flov?- along 
the surface* The pressure is computed directly from the velocity 
at the points of the ellipsoid, now supposed -ic be stationary in 
the flowing fluid. For this is a steady flow, an^i hence Bernouilli'e' 
equation for the pressure holds true, vXa.' p + P = corst. 
That is, the pressure is equal to an arbitrary constant pressure 

2 P 

minus where V^^ denotes the velocity. The points of 

greatest velocity are those of smallest pressure or of greatest 

suction* The curves of eqtial velocity are also the curves of equal 
» pressure. 

In practice, we are chiefly interested in rather elongated 
ellipsoids of rotation, and the angle ot, between the principal _ 
axis and the direction of motion is small* W.ith very elongated 
ellipsoids, ki is about 1 and ka is very small. Hence A is 
about 3 and B is about 1, and the angle between the direction 
of the line of maximm velocity and the axis is about twice as 
large as the angle between the direction of motion and .the axis. _ 
The maximum velocity is always greater than the velocity of motion. 

The difference between the largest negative pressure and the pres- 

sure in the undisturbed atmosphere is 

V^|-^ (l + cos= a + (1 + kgf sin^ a - 1 | 



